INFINITE MATRIX REPRESENTATIONS OF ISOTROPIC 
PSEUDODIFFERENTIAL OPERATORS 



OTIS CHODOSH 

O 

Abstract. We characterize isotropic pseudodifferential operators (elements 
$-H 1 of the Shubin calculus) by their action on Hermite functions. We show that a 

£L[ continuous linear operator A : <S(R) — > S'(R) is an isotropic pseudodifferential 

operator of order r if and only if its "matrix" (ffW) mj0 := (A<f> n , cj>m) L 2 ( R j is 
rapidly decreasing away from the diagonal {m = n}, order r/2 in m + n, where 
applying the discrete difference operator along the diagonal decreases the order 
by one. Here <j> m is the m-th Hermite function. As an application, we give an 
isotropic version of the Beals commutator characterization of pseudodifferential 
operators, showing that if we define H = —d^+x 2 to be the harmonic oscillator 
and Z the map extended linearly from Z{<f>k) = 4>k-i, then a continuous linear 
operator A : S(K) — ¥ <S'(R) is an isotropic pseudodifferential operator of order 
r if and only if commuting A a times with H and fi times with Z results in an 
bounded linear operator H?+ r ~ 2p (M.) -> Hf (R), for all s e R and a,/3 S N . 
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1. Introduction 

In this paper we examine the action of isotropic pseudodifferential operators on 
the Hermite functions. We give a necessary and sufficient condition for an operator 
to be an isotropic pseudodifferential operator based on the form of action. This 
result is a consolidation of results originating in the author's honors thesis, [2], 
which contains, in addition to the material given here, a proof of a similar theorem 
for pseudodifferential operators on a torus (which seems to be well known, although 
a proof could not be located in the literature) as well as more detailed background 
about Hermite functions and isotropic pseudodifferential operators. In addition to 
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the above result, we give an application of our discretization, proving a commutator 
characterization of isotropic pseudodifferential operators in the sense of Beals. 
In order to state our main theorem we first give the following definition: 

Definition 1.1. If we define the discrete difference operator A on a function K : 
N x N -> R by 

(AK)(m, n) = K(m + 1, n + 1) - K(m, n) 

(writing A a to signify applying the difference operator a times), then we will say 
that a function K : No x No — > R is a symbol matrix of order r if for all a, N £ No, 
there is C Qj at > such that 

\{A a K){m,n)\ < C a , N (l+m + n) r - a (l + \m-n\)- N . 

We will denote the set of symbol matrices of order r, SM r (No). 

With this definition, we can state our main result: 

Theorem 1.2. A linear operator, A : <S(R) — > <S'(R) is an isotropic pseudodiffer- 
ential operator of order r (as defined in Section [3|) if and only if the "matrix of 
A" 

K {A) : No x N -> R 

defined by (m,n) H ► (A(f) n , 4>m)L 2 (M.) * s an order r/2 symbol matrix (where <j) n is the 
n-th Hermite function, as defined in Section^. 

Perhaps the simplest example of this theorem is when A is a power of the har- 
monic oscillator. It can be shown that for s e R, (1 + H) s / 2 e *f so (R) is an elliptic 
operator, where H is the harmonic oscillator, H = +x 2 (see, for example, [9], 
Theorems 11.10. 1 and II. 11. 2 concerning powers of elliptic operators). Furthermore, 
because {\+H)<j> n = (2 + 2n)<p n , it can be shown that (l + H)^ 2 c/) n = (2 + 2n) s / 2 0„. 
Thus 

(i^ 1+ff ) s/2 ) m ,n = ((1 + H) s ' 2 <p n , <f> m ) L2 = (2 + 2n) s / 2 S m , n . 

Here we have written {K^ 1+H ^ 1 ) m ,n when we really mean K^ 1+H ^ 1 (m,n), a 
shorthand we will use frequently below. Notice that this is a symbol matrix of 
order s/2: because (K^ 1+H ^ 1 ) m . n = for m ^ n, it is enough to show that 

\ {A a K a+Hy /2 )nn \ < Ca (i + n) s/2 - a . 

This follows from the fact that d^(l + x) s < (l + a;) s_fc and the mean value theorem. 
Thus, powers of the harmonic oscillator, which are some of the simplest examples 
of isotropic pseudodifferential operators satisfy the statement of the theorem. 

Theorem 11.21 is similar to the situation for pseudodifferential operators on the 
d-dimensional torus T d . In this case, it is well known that for an operator A : 
C°°(T d ) — > C°°(T d ), A is an order r pseudodifferential operator on T d if and only 
j£K^ A '(m,n) := (A(e m ' x ),e lm ' x ))L 2 is a symbol matrix of order r (where we extend 
the definition of symbol matrix to a map K : Z d x Z d — > R in the obvious manner, 
letting the difference operator act in any of the d-directions). For the dimension 
d = 1, this is very similar to the isotropic case discussed in this paper. However, 
for higher dimensions, there is not such a simple characterization of the matrices of 
^•" so (R d ). We discuss this further in the end of Section HI A proof of the theorem 
on the torus can be found in [2]. 
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Related forms of discretization of pseudodifferential operators have been studied 
by various authors. Ruzhansky and Turunen in [6j|8] have considered Fourier series 
discretization of toroidal pseudodifferential operators . In their work, they consider 
"quantized" symbols of the form a(x, £) £ C°°(T d x Z d ) with 

\d^A^a(x,O\<C a Al + m m - 

where A' 3 is the hnite difference operator applied (3 times in the £ variable. These 
symbols then define a toroidal pseudodifferential operator A a £ ^i m (T n ) via the 
action on test functions 4> £ C°°(T d ) given by 

A a 4>(x)= V / e< x -y^a( Xl ^{y)dy. 

They show that any pseudodifferential operator on the torus has such a symbol, 
and describe how the quantized symbols relate to the usual notion. They go on 
to discuss Fourier intergral operators on the torus and applications to hyperbolic 
equations. Their quantization is somewhat different than the one described in this 
paper, because they are quantizing only in the £ variable, and not in the x variable 
as well. In addition, we remark that in [7] , they have extended this quantization to 
pseudodifferential operators over general compact Lie groups, and it is interesting 
to note that they have also observed a different form of off diagonal decay of the 
symbols in the setting of e.g. matrix valued symbols for operators in \p m (SU(2)) 
(cf. Chapter 12). 

Finally, our characterization of isotropic pseudodifferential operators allows us 
to give the following commutator characterization in the sense of Beals [l . We 
define the operator Z : 6>(R) — > <S(R) by 4>k i-> 4>k-i (extending linearly), and 
recall that the harmonic oscillator is H = —■§-* + x 2 . and from this, we define the 
commutation operators Z(A) := [A, Z) and H(A) := [A, H], allowing us to state 

Theorem 1.3. An operator A : S(R) -> S'(R) has A £ *- so (K) if and only if for 
all a, P £ N and s £ ffi 

(i.i) h^{z^{a)) £ L{H r + s - 2 ' 3 m -> Hum)- 

This theorem gives an isotropic version of the classical theorem of Beals, which 
says that if commutators of a linear operator with Xj and dj have appropriate map- 
ping properties, then it is a (standard) pseudodifferential operator. It is important 
to note that the above operator Z is an isotropic pseudodifferential operator, but 
not a differential operator. 

In the following, we fix a choice of Fourier transform, defining 

/•./; O = /(0 = / e~**f{x) dx. 
With this convention, the inverse Fourier transform is 

r-'m = I e^f(o 

where d£ := {2i:)- d d^ 

Additionally, we will also use the standard notation f(x) < g(x) to denote that 
there is some C > such that f(x) < Cg(x) for all x. When both sides of the 
equation depend on multiple variables it should be clear from context which ones 
C depends on. We also write <S(R) for Schwartz functions (/ : R — > R such that 
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for all a, /3 G N = {0, 1,2,...}, we have that ||x Q ^/|| L 2 (R) < oo). We further let 
<S'(R) be the space of tempered distributions, the topological dual to the Schwartz 
space <S(R). 

We have compiled various elementary results about symbol matrices in Appendix 
IA1 most of which are used in the proof of our main result, Theorem ll.2l Appendixes 
iBl and [C] contain proofs of Lemmas which are needed in the proof of the main 
theorem, but we felt were somewhat extraneous to the exposition, as they both 
are somewhat messy calculations. Appendix [C] could potentially be of independent 
interest, in it we prove that for A G ^^(R) (for eq < 1/2), if A is bounded 
L 2 — > L 2 , then the symbol of A is bounded. As such, it is a partial converse to the 
L 2 boundedness of order pseudodifferential operators. 



Much of this paper represents a part of my undergraduate honors thesis. I am 
very grateful to my advisor, Andras Vasy, for teaching me about pseudodifferential 
operators and microlocal analysis, and for suggesting and then greatly assisting me 
with this work. I would like to thank Richard Melrose, Michael Ruzhansky, and 
Michael Taylor for their helpful comments and suggestions. I would also like to 
thank the referees for their careful reading and comments. Part of this work was 
completed while supported by a Stanford University VPUE research grant. 

3. Isotropic Pseudodifferential Operators and Hermite Functions 

Isotropic pseudodifferential operators, also known as the "Shubin class" in the 
literature, are defined to be pseudodifferential operators with symbols a(x, £) G 
C°°(R 2 ) such that for all a, /3 G N , there is a C ai p > so that for x, £ G R 



Here we follow Melrose's development of the class (c.f. [4]). Note that these differ 
from "regular" symbols in two ways, namely that we require decay in both x and 
£ (regular symbols only require decay in £) and we also require x derivatives to 
improve the decay (for regular symbols the x derivatives are only required to not 
make the decay worse). We denote the above class of symbols by S^ (R), and class 
of operators with these symbols \&f so (R), which we call isotropic pseudodifferential 
operators of order r. 

We recall that given an (isotropic) symbol a(x, £), we obtain an (isotropic) pseu- 
dodifferential operator A by setting 



for Schwartz functions u. In, e.g. [4], it is shown that this is a well defined continuous 
operator on Schwartz functions. Additionally, it can be shown that the composition 
of two isotropic pseudodifferential operators is again an isotropic pseudodifferential 
operator, and that the adjoint of an isotropic pseudo-differential operator is again 
an isotropic pseudo-differential operator of the same order (see, e.g. Melrose [3] 
Proposition 4.1 and Theorem 4.1). Furthermore, it can be shown that elliptic 
elements of $r (R") have two sided parametrices in \J/r£(K d ). Additionally, because 
^•^(R") C ^°(R ), we have L 2 boundedness of order isotropic pseudo-differential 
operators. 
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We will need the notion of an elliptic isotropic pseudodiffcrential operator of 
order r, which is is an isotropic pseudodifferential operator whose symbol obeys 

\a(x,£)\>C(l + \x\ + \£\y 

for x, £ outside of some compact set. We further remark that, as in the case of 
standard symbols, the full symbol is not uniquely determined by the operator, but 
its image in S^ /S^~ is uniquely determined, and we call this image the principal 
symbol of the operator A. 

One property of isotropic calculus that is different from the standard calculus 
is closure under the Fourier transform in the following sense: defining A by the 

formula Au = Au, it can be shown that A is a isotropic pseudodifferential operator 
of the same order with principal symbol a(x, £) := a(£, —x). 

Finally we discuss the Hermite functions on M. d . If we define the j-th creation 
operator Cj = (— g|- +Xj) , and let </> = -^rfze~\ x \ / 2 , then to define the higher 
Hermite functions (depending on a multi- index a), we inductively let 4> a +ej — 

1 = C]4> a . These are an orthonormal basis for L 2 (M. d ) and also span a dense 

subset of S(R d ). Note that H<j) a = (1 + 2\a\)(j> a (where H = - Y? j=1 J^r + \x\ 2 

is the harmonic oscillator), and J 7 </>&(£) = (— i) Q (27r) d / 2 0fc(£), so in particular the 
Hermite functions are eigenfunctions of both the harmonic oscillator and the Fourier 
transform. More properties and proofs of these claims can be found in many places, 
e.g. Chapter 6 of [10] , Alternatively, the results that we will make use of are 
described in detail in Chapter 4 of [2]. For later use, we record the following 
(c.f. ,2} Lemma 4.8]) 

Lemma 3.1. For a multi-index n £ Z d with > for 1 < k < d we have 



(3-2) x k (p n {x) = y Uk ^ - <t> n+ek (x) + J Y<t>n-e h (%)■ 

4. Proof of Theorem 11.21 

In this section we prove Theorem 11.21 showing that order r isotropic pseudodif- 
ferential operators are the same as operators whose associated matrices are symbol 
matrices of order r/2. 

Proof that A e *^ S0 (M) => e SM r / 2 (N )- We assume that A G 1Tso( R ) is S iven - 

We will use L? boundcdncss of order pseudodifferential operators to show that 
the matrix of A is of order r/2 in m and n. Then, an integration by parts argument 
will show that the matrix is rapidly decreasing off of the diagonal. Finally, we will 
show that applying the difference operator gives a matrix that is comparable to the 
matrix of an operator one order lower, which will complete the proof. 

It can be shown (as discussed in the introduction) that if H = — + x 2 is the 
harmonic oscillator then (1 + H) s / 2 € ^^(K) with principal symbol (1 + \x\ 2 + 
|£| 2 ) s / 2 . For a proof of this, see gj, Theorems II.10.1 and II.11.2. We define 

(4.1) A t ={l + H)~ rt ' 2 A{l + #)-Ki-*)/2, 
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Because A t G *° so (E), it is bounded on L 2 (R). This implies that 

\(K {At) ) m , n \ = \(A t (f> n ,<f> m )\ < \\A t \\ c{L 2 m < 00. 

By Lemma IA.1| the matrix of a product of two operators is the product of the 
matrices of the operators. Thus, using the convention that repeated indices are 
summed, 

{K^ ] ) m . n = ((i + H)- rt / 2 ^> m ) L2(T , )) (^)) feJ ((i + ffr r(1 ~ t)/2 ^,<M^(T*)) 

= £ fc ,m(2 + 2n)- rt / 2 {K^) K] 5 h „(2 + 2m)- r ' 1 - i »/ 2 
= (2 + 2n)- rt / 2 (2 + 2m)- r ^-^ 2 (K^) mtn 
= C(l + n)-^ 2 (l + m)-^ 1 -*)/ 2 ^^))^^ 

Combining this with the above bounds on (K^ At ') min (because it is the matrix of 
a order pseudodifferential operator), we thus have that 

If the order of A, r, is positive, taking 4=1/2 gives 

i(^ (A) ) m ,«i<(i+mr /4 (i+n) r/4 

= (1 + m + n + mn) r / 4 

< (1 + 2m + 2n + 2mn) r/4 

< (1 + 2m + 2n + 2mn + m 2 + n 2 ) r/A 
= (l + m + n) r / 2 

and if the order r < 0, taking t — 0, 1 gives 

|(^ (A) ) m ,n| < (l + min(m,«)) r / 2 

< (l + 2- 1 (m + n) r/2 

< (l + m + n) r / 2 . 

In order to show off diagonal decay of K^ A \ we first show how to simplify 
the integral expression of (K^) m , n by using the fact that Hermite functions are 
eigenfunctions of the Fourier transform: 

[KW) m , n = {A<j> n ,ct> m ) 

e l ( x -yKa(x,OMy)<f>m(x)dyazdx 

e lx( -a{x, QFy^t [<t>n] {Q<t>m (x)d^dx 
e l ^a(x,C)[V2^(-i) n M0<t> m {x)d^dx 
M-i) n I e^a(x, 0<Pn(0<t>m (x)d^dx. 



Letting H x = x 2 — J^p and = £ 2 — , we thus have that 
= 2(m+l-n-l)(K^) m , n 
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= v^F(-i) n / e lx <a(x,0(H x - JT £ ) [</»„ (£)<M*)] dxd£ 

JR d xE d 

- V^F(-i)" / (H x - H ( ) [e ix <a(x, 0] M04>M dxd^. 

JR d xM d 

Now, we will show that the difference of the harmonic oscillators applied to e lx ^a(x, £) 
results in e lx ^a(x, £) where a(x,£) is a new symbol, still of order r. 



(H x -H s )[e^a(x,0] 



d 2 

d 2 x <9 2 £ 



(x 2 -e+e- ^ " ixi 



+ e 



e*°*a{x,Z) 
a(x,£) 

d d 2 



^ dx ^ X d£ d 2 x dt; 2 



e 

= e^a(x,0- 

Applying this repeatedly, we see that for any N > 0, multiplying K^n by (m — n) N 
is the same as considering (K^) m ^ n for some A e ^[ so (E), and thus, repeating 
the boundedness argument given above for A, we thus have 

(4.2) \(K (A) Un\ < CW(1 + \n- m\)- N (l + m + nf' 2 . 

It remains to be shown that the difference operator lowers the degree of the 
matrix. For the difference operator, notice that (letting C := -^(x + d x ) be the 

annihilation operator and = -^{x — d x ) the creation operator) 

(Aif (A) ) TO; „ =(^„ + i,<^ m+ i) - (A(p n ,(f) m ) 

-.{AC^ n ,CU m )-{A4> n ,(t) m ) 



d 



d d 2 



d 2 



' 2< dx + 2ix dt; d 2 x + d 2 ^ 



a(x,£) 



y/(m + l)(n+l) 
1 

V(m + l)(n+TJ 
1 

^(772+ l)(n+ 1) 

/ n + 1 



y / (m + l)(n+ 1) 
1 

+ v /(m + l)(n = +iy 



(CAC*^,^) - (A(j> n ,(j> m ) 

((AC+[C,A])CUn,<f>m) - (A<f> n ,<f> m ) 

1 (A(f> n ,<p m ) 



(([C,A])CUn,<f>m) 



n+l 
m + 1 



1 (A(j> n ,<j> m ) 



x /(m + l)(n + l) 



(([C,A])Ct0 n ,0 m ) 



Vri+T — vm + 1 
Vm+ 1 



{A(j> n ,(t> m ) 
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+ = (([(?, A])C^„, ^ 

a/ (to + l)(n + 1) 



n — to 



V / roTT(v / nTT + Vm + 1) 
1 



(A/>„, 0„ 



V(m + l)(n + l) 



(([C,A])C^ n ,^ 



Because 



(\/m + 1 + y/n + l) 2 = to + n + 2 + 2^(m + l)(n + 1) > (1 + \m\ + \n\) 

and 

(to + l)(n + 1) = (1 + to + ?t, + mn) > (1 + to + n) 2 , 

combining these inequalities with (|4.2[) (the boundedness argument applies for 
[C, A]C\ which is a degree r + 1 isotropic operator) we see that the difference 
operator will lower the degree in to + n by 1 and still preserve the off diagonal 
boundedness. Repeating this argument inductively for higher powers of the differ- 
ence operator, we see that K^ A > is a symbol matrix of order r/2. 

□ 

Proof that £ SM r / 2 (N ) A <= *- so (M). We assume that A : 5(R) S'(M) 

is given, with if^) S SM r ' 2 (No)- By multiplying by appropriate powers of the 
harmonic oscillator, we will show that we get a Hilbert-Schmidt operator. Using 
this, we will show that we have weighted L 2 bounds on the derivatives of the symbol 
of A. This will give weighted L°° bounds on derivatives of the symbol, which will 
show that A g ^\Z, e (R) for all e > 0. Finally a sort of converse to L 2 boundedness 
of order operators will show that we can take e = 0, as desired. 
Let K(x, y) be the Schwartz kernel of A. We define a distribution 



This definition gives that (where in the sequel we will follow the standard practice 
of writing an integral where we actually mean pairing with a distribution) 

However, if we replace a(x, £) with £) we see from the following calculation 

that this results in a symbol matrix of order (r — l)/2 (we will use Lemma 13. f I 
which gives a formula for derivatives of Hermite functions as well as the product of 
Hcrmite functions with a linear function): 

M-i) n f e^^(x,0M0^m(x)d^dx 
Jrxm ox 

/ 2^(-i) n+1 [ (t + i-p) [e lx <a(x,0] MO0m(x)dlidx 
Jlxl V ox J 

^(-i) n+1 [ e ix -ta(x, U - MO0m(x)dZdx 
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-(t>n-l{d4>m{x) +ij <p n (^)(j) m +l(x) 

-ij -j(t>n(Q<f>m-i (x) J d£dx 



m,n+l \/ 2 m ' n— 1 V 2 w+lj^* y 2 m ~ li n 



( A) / n + l ( A) /n+1 (A) 

"-"■m.n+l \/ o n m-l,n 'V 2 m— l,n 



It is clear that the order of symbol matrices is multiplicative and that the difference 
operator applied to a symbol matrix is a symbol matrix of one order lower, which 
shows that the terms with the difference operator are symbol matrices of order 
(r — l)/2. Lemma I A. 2 1 (which deals with multiplying a symbol matrix by a func- 
tion which is not necessarily decaying off the diagonal, but does have appropriate 
diagonal decay) shows that the remaining terms are also symbol matrices of order 
(r — l)/2, so thhe sum of all four terms is also a symbol matrix of order (r — l)/2. 
A similar calculation holds for £ derivatives. 

We let (for a, (i € No) A a '^ be the operator with Schwartz kernel 

From our above results, we have A a '$ : S(R) -> S'(R) and its matrix is a symbol 
matrix of order (r-a- (3)/2. Now, for s e R, we define Bf$ = (1 + H)- s / 2 A a ^. 
By Lemma fA.31 (which says that order r < —1/2 symbol matrices are £ 2 summable, 
implying that the corresponding operator has an L 2 summable kernel) if r — s — 
a — (3 < — 1, then B™'* is Hilbert-Schmidt because its kernel is square summable 
(see e.g. Theorem VI. 23 in Reed-Simon [5]). Thus, for the kernel of B"^ we have 
K?'P(x,y) e L 2 (R x R). This implies that for r(x,g) the symbol of (1 + H) s / 2 , 

A a ^f = (1 + Hf' 2 ( f K^(y, z)f(z) dz 

\JR J 

e^- y K(x, (V, *)/(«) dzdy<%. 



Taking the Fourier transform 

d a+(s a 



() ,,^(^)e(i + w 2 + iei 2 r /2 i 2 (Rx 
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for r — s — a — < — 1, because the Fourier transform of Kf'P is in L 2 (R x R), 
and the principal symbol of (1 + H)~ s / 2 is (1 + x 2 + £ 2 )~ s / 2 and the lower order 
terms will only improve the integrability. By Sobolev embedding, this implies that 
actually a{x,£) G C*°°(R). 

We have weighted L 2 (R x K) bounds on the derivatives of a, but to show that 
a is a symbol, we would like weighted L°°(R x R) bounds on the derivatives. This 
is a standard argument, we include it in Appendix [B] More precisely, Lemma IB. II 
shows that for r' with r — r' — a — /3 < 



d a+/3 a 
dx a d^ 



(i,0 g (i + N 2 + |e| 2 ) r ' /2 r(ixt). 



Thus, we see that A G $f s ' (R) for r' > r, and B^* 3 G ^^ s ~ a ~ (SL) for r' > r. 
Thus, to complete the proof, all that we must show is that A G 4'-l (R). Notice 
that this implies that for r — s = a + (3, B^' 13 G ^^(K) for any e > 0. However, 
since the symbol matrix of -B" ,/3 is order 0, by Lemma IA.5I B*'P is bounded on 
L 2 (R). Thus, we see that for r - s = a + /3, is bounded on L 2 (R). Using the 
following lemma, we can thus conclude that the symbol of Bf'° is bounded, which 
then implies that A £ 5 , f so (R) because that will give the necessary bounds on each 
derivative. 

Lemma 4.1. For Q G *-° (R) for some e < 1/2, if Q is bounded on L 2 (R), then 
the symbol of Q, ctl(Q) = q{%,£) is bounded, i.e. 

sup \q(x,£)\ < oo 

The proof of this is a calculation based on the idea of letting a pseudodifferential 
operator act on an exponent to gain control of the symbol. We give it in Appendix 

o □ 

This theorem does not hold as stated in R d for d > 2. For example, consider C\, 
the creation operator in the first coordinate. Recall that 

C\4>n = y/fli + l^n+ei 

Thus 

(X (Clt) ) m ,„ = {C\(j) n ,4>m) = \fni + lS n+eiim 

This has the appropriate decay; away from the diagonal it is zero, and for n+ei = m 
it is order 1/2 in as in the d — 1 case. However, taking the difference operator 
in the first coordinate gives 

2 , 



y/m + 3 + y/m + 1 

This is order —1/2 in ni, but merely bounded for e.g. n^, so in general the order in 
\n\ is only 0. Furthermore, taking further difference operators in the first coordinate 
do not change the order in |n| at all. 

It is possible to prove an analogous theorem for dimension d > 2 by using 
arguments similar to the above proof. 
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Theorem 4.2. For A a continuous linear operator S(R d ) -> S'(R d ), A € ^ (R d ) 
if and only if (K^ A ') mi „ has the following property for all multi-indexes a, (3 G Nq 

(4.3) \{U a ^ K {A) ) m , n \ < (1 + \m\ + \n\)( r -^~WV 2 

where we define (D Xk K) nrn by 



n k + l 



K 



n-\-ek,m 



and (D (k K) nm by 



rrik + 1 



K, 



Kn 



m k + 1 



Kn 



n k + l 



K„ 



K„ 



id let = (D xi r ■ ■ ■ (a Xd r ny 1 ■ ■ ■ 

The n Q ^ operator is really the matrix analogue of applying d x d^ to the symbol 



of A, as we illustrate in the case of (a,/3) = (efc,0): 

nm 



n k 



1 



-K n 



K n . 



m k + 1 



K n , 



m+e k 



K 



= {2n) d / 2 {-i) n+1 



a{x,£) 



n k + 1 



\n+l 



2 

{2ix) d / 2 { 
{2n) d / 2 {-i) n+1 
(2n) d/2 (-i) n+2 



m k + 1 



„ d£dx 



dxi 



4>n(£,)(t>m{x)d£,dx 



d 



^ k + i dx~ k ~ ' e ix *a(x,Z)<j> n (£)(j) m (x)d£dx 



da 

dx k 



{x,£,)<t>n{£,)<i>m{x)d£,dx 



The proof of Theorem 14.21 follows the one dimensional case almost identically, 
except that the above calculation replaces the difference operator calculations. For 
example, to show that for A £ \&r o (R d ), equation (|4.3p is satisfied, the argument 
for (a, (3) = follows in exactly the same way from L 2 boundedness of zeroth or- 
der isotropic pseudodifferential operators. Then, applying the box operator gives 
a matrix which is the matrix of an isotropic pseudodifferential operator with the 
appropriate derivatives on the symbol, so again we have the desired bounds. Con- 
versely, it is possible to show (as we discuss below) that matrices obeying (|4.3[) with 
?' < —d/2 are square summable, and thus correspond to Hilbert-Schmidt operators. 
Thus, we can establish weighted bounds on the distribution that we would like to 
show is the symbol (exactly as in the one dimensional case). Finally, derivatives 
of the "symbol" correspond to applying O a '@ to the matrix, and we thus have 
weighted bounds on the derivatives of the symbol, allowing us to conclude as in one 
dimension. 

The conditions of (|4. 3[) are somewhat unsatisfying, as they are considerably more 
complex than the simple difference operator results in the one dimensional case. For 
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example, it can be shown that infinite matrices obeying (|4.3|) are rapidly decreasing 
off of the diagonal, which is not obvious from the condition. To see this, consider 
the following identity 

'm k Ink 

— {u Xk + u & )K n _ m _ ei! - w -y' 



(m fc - n k )K nm = J — (D Xk + U ih )K n ^ m - ek ~ J — (D Xh - U^ h )K n - eh , m . 



It is not hard to see that the right hand side still obeys (|4.3p with the same r, which 
shows that K nrn is rapidly decaying off of the diagonal (in particular, this shows 
that condition (|4.3p with r < —d/2 implies that the matrix is square summable, as 
promised above). 

It is likely that the issues relating to the difference operator in higher dimensions 
is a consequence of the high multiplicity of the harmonic oscillator eigenspaces. It 
seems possible that a better result could be obtained through a "rearranging" of 
the eigenspaces. That is, the difference operator compares the operator's action on 
4> n and <f>n+e k j elements of the \n\ +d and the \n\ + d+l eigenspaces, but there is no 
reason that these are the proper elements from these two eigenspaces to compare. 
However, we have been unable to find a satisfactory manner in which to compare 
the two eigenspaces in order to obtain a simpler condition on the matrix. 



5. An Isotropic Beals Theorem 

We begin by recalling the classical theorem of Beals characterizing (standard) 
pseudodifferential operators on R d . For an operator A : S(R d ) — > 6>'(R d ), we 
define Lk(A) :— [xk,A] (where [-B, C] := BC — CB is the commutator) as well as 
R k {A) := [dig, A}. Theorem 1.4 in [I] gives 

Theorem 5.1 (Beals commutator characterization). A continuous linear map 

A : S{R d ) -> S'(R n ) 

is a pseudodifferential operator of order r, i.e. A € ^ r (R d ) if and only if for a, (3 £ 
(No) d we have that 

LW(RW(A)) £ £(H r+ W-W(R d ) -> L 2 (R d )). 

Here, we have written C(H r+ ^~^(R d ) -> L 2 (R d )) for bounded linear maps from 
the r + |/3| — |a| Sobolev space to L 2 and have written for L composed with 
itself a times (similarly for R^' ). 

By using Theorem II. 2 [ we can give a version of commutator characterization in 
the setting of isotropic pseudodifferential operators. First, we recall the definition 
of isotropic Sobolev spaces 

tffsoW = {/ € S'{R) : (1 + H) s ' 2 f £ L 2 (R)} 

= lf£ S'(R) : J2(l + k) s (/, 4>kf L , m < oo 1 
I fc=o J 
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We define the operator Z : S(M.) — > <S(R) by fa H> fa~\ (extending linearly). It 
is clear that Z^(fa) — fa+i. Furthermore, notice that the matrix of Z is 

/0 \ 
1 

1 

1 



K (z) 



\ 



J 



and thus, by Theorem 1 1.21 Z G ^P so (R) (but it is certainly not a differential opera- 



tor). It is clear that ZZ> 
has matrix 



Id. Also, recall the harmonic oscillator H = —d% 

( l \ 

3 

K (H) 



Given these two operators, we define Z(A) := [A, Z] and H(A) := [A, H\. Now, we 
recall our version of the isotropic Beals theorem and give the proof: 



Theorem 11.31 An operator A : S{ 
all Q,/3gN and sgl 



S' 



with Aef! 



if and only if for 



(5.1) 



Proof. It is enough to show that (|5.1[) for all a, /3 e N implies that A G *[ so (R). 
We will show that if ( A > is a symbol matrix of the appropriate order, G 
SW/ 2 (N ), and then use Theorem II .21 to conclude the desired result. 

First, notice that by repeating the argument following equation (|4.1[) (conju- 
gating by the appropriate powers of the harmonic oscillator) we have that if an 
operator B : S(R) -> «S'(M) with B G C(H[+ S (R) -> #f so (R)) for all s6l and 
some r G M, then 

l(# (B) )m,„| <(l + m + nf/ 2 . 

In particular, taking B = A in the above, shows that if has the desired overall 
decay of order r/2. It thus remains to show that the finite difference operator 
reduces this order of decay and off diagonal decay, following assumption (|5.1|) . For 
the finite difference operator, notice that 



(Aif (A) ) mj „ = (A0„ + i,0 m+ i) L2(]p 

= {AZU n ,ZUrn) L2{ 

= (ZAZ^ fa,fa % ) L2{1 
= ((AZ-[A, Z])Z^fa,cp m ) L2 
= ({A-[A,Z]Z^)fa,ct> m ) LHM 

= - ([A Z \ Z] faAm) L 2 (m) ■ 



(Afa,(f) m ) L 2( R ) 

- (Afa,(f) m ) L 2( R ) 

(A fa, <fi m ) L 2( R ) 

- (A0„,0 m ) L2( 

(A fa, 0m)i2( R ) 



Because K ^ ) G 5M°(N ) is invertible, modulo a finite dimensional error (because 
Ztf = Id but Z^Z{fa) = fa for k > and for k = 0), dSUJ) and the argument 
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above with B = Z(A) (along with the multiplicative property of kernel matrices, 
Lemma I A. 1|) implies that 

| (A ^U)) m>n | <(i + m + Tl f/2-i. 
Repeating the above shows that 

|(A" K^) m , n \<(l+m + nyl 2 -P. 
For the off diagonal decay, notice that for any B £ £(H? so (M) —> L 2 (M)) 
2(m - n){K ^) m ,„ = (1 + 2m - 1 - 2n)(K {B \ n , n 

= (B(f> n , (1 + 2m)(j) m ) L2(R) - (5(1 + 2n)<t> n ,<t>m) L 2 (R) 
= (B4> n ,H<f> m } L2 ^ - (BH(f) n ,<f> m ) L 2( M ) 

= (HB<f> n ,(() m } L2 g^ - (BH(f) n ,<f> m } L2 Q^ 

= - ([B,H](f> n ,(f) m ) i2 ( R ) • 

Thus, repeating this inductively, and combining with the above argument, we see 
that for all N > 

(m - n) N \(A? K^) m J < (1 + m + nfl 2 ^ 
which allows us to conclude K ^ € SM r / 2 (N ), as desired. □ 

Appendix A. Symbol Matrix Properties 

In this appendix we compile various properties about matrices of operators and 
symbol matrices that were used in the preceding sections. First, we show that 
associating a map with its matrix is multiplicative in the same sense as finite di- 
mensional linear operators. 

Lemma A.l. The map from operators A : S(R) — > <S(R) to matrices K^ A ' : 
No x No — > K is a homomorphism where we define the product of two matrices to be 
the infinite dimensional analogue of matrix multiplication. Namely, for continuous 
linear A, B : S(R) -> S(R) 



Proof. For any Schwartz function / € 5(R), we have that the following sum con- 
verges in 5(R): 

^2 , ? !) fc(/'^)i 2 (R) / 

0<fc<JV 

as N — > co. Now, using this fact, along with the fact that A, B are continuous as 
operators on 6>(K), the following calculation is justified: 

(K (AB) ) m ,„ = (AB(j) n , 4> m ) L2(R s 



\ VfcGNo / 

^(K^k^A^),^ 
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E {{K^) k ,nA{4> k U;, l/Liii 



□ 

Now we show that symbol matrices are closed under multiplying by a function 
which is polynomially bounded with the difference operator lowering the decay by 
one order. 

Lemma A. 2. If f(m,n) is a symbol matrix of order r and g(m,n) is any function 
g : No x No -> K satisfying 

|(AT ff ) (m,n)\ < C 7 (l + \m\ + \n\) s -^ 

for 7 G N and seR, £/ien /(to, n)g(m, n) is a symbol matrix of order r + s. 

Proof. It is easy to see that from the conditions on /(m, n) and g{m, n) we have 

\[fg](m, n)\ < C' Nfl {\ + \n- m\y N (l + \n\ + \m\) r+s , 

so the only thing that is left to check is that applying the difference operator gives 
improved decay along the diagonal, which follows from 

A [fg] (to, n) = A [/] (to, n) g(m + 1, n + 1) + f{m, n) A [g] (to, n). 

Repeating the above calculation finishes the proof. □ 

Because of the rapid decay of a symbol matrix off of the diagonal, we expect 
many of the properties of a diagonal matrix to hold for a symbol matrix. One such 
property is € 2 (No x No) membership for symbol matrices: 

Lemma A. 3. If r < —1/2, and K £ 5M r (N ) is any symbol matrix of order r 
then {{K) m>n } € £ 2 (N Q x N ). 

Proof. For a fixed < s < 1 and N e No, we bound the sum as follows 

m,n£No neN .fe> — n 

= J2 \(K) n+k ,n\ 2 + E \( K )n+k,n\ 2 

|fe|<|n| s ,fc>-n |fe|>|n| 5 ,fe>-rt 

< (l + n + n + kf r 

|fc|<|ri| s ,fc>-n 

+ E (i + ifcir 2JV (i+™+™ + fc) 2r 

|fe|>|n| s 

< J2 a+-) 2r + E (i + ifcir 2Ar (i + ifci) 2M/s 

|fe|<|ri| s ,fe>-n |fc|>|n| s ,fc>-n 

< E #{ fc: 1^1 < \n\ s }(l + n) 2r 

n€N 
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+ E / (i + Mr 2W+2r/s dx 



<EH s (i + «) 

n£N 



2r 



/•oo 

+ E / (l + t)" 2W+2r/s dt 

E H s (i + ^) 2 



< \ . . - / ( i ^ ^ 2 / 

n£N 

+ E( 1 + M s ) 1 ~ 2W+2r/s - 

neN 

For s + 2r < — 1 and s(l — 2iV) + 2r < — 1 the final terms will be finite. Because 
r < —1/2, we can choose s > such that 2r < — 1 — s. For this s, we can then 
choose an N large enough such that s(l — 27V) < — 1 — 2r. This completes the 
proof. □ 

Given an arbitrary symbol matrix, it defines an operator on Schwartz functions, 
as is shown in the next lemma. 

Lemma A. 4. If f(m, n) is a symbol matrix of order r, then it defines a continuous 
operator 5(E) -> 5(E). 

Proof. It is clear that our desired result is equivalent to proving that the operator 

F : {x fe } feeNo i-> E f(n,k)x k 
fceN 

is a continuous linear operator on the space of rapidly decreasing sequences 

«(N ) - {{x fc } feeNo : |(1 + \k\) N x k \ < C N VN e N } 
where the topology on s(No) is given by the seminorms 

\\xk\\N = sup |(1 + k) N x k \. 

k 

Because a symbol matrix of negative order is also a symbol matrix of zero order, 
we may assume that r > 0. Notice that multiplication by (1 + \k\) s is a continuous 
linear bijection s(No) — > s(No). Denote this operator by L s . Choosing any N, 
notice that 



\L S o Fo Li({x k })\\ N 

J2(l + \n\) s f(n,k)(l + \k\) l x k 



fceNo 



N 



sup(l + \n\) 



N+s 



E f(n,k)(l + \k\y 



■l- k 



< 



< 



Csup(l + n) N+s E (1 + M + |fc|) r (l + |fc|)'|a*| 



fceNo 



C||x fc || sup(l + M)^ E (1 + H + WH 1 + 
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< C\\x k \\ sup(l + \n\) N+ ' J2 C 1 + MH 1 + IW 1 + l fc D' 

" fcGN 

= C\\x k \\ sup(l + |«|)^+-+- £ (1 + 

" fcGN 

In particular, taking I < — 1 — r, and s < —N — r, we have that 

\\L a oFoLi({x k })\\ N < Ci, a \\x k \\ . 
Thus, for these s, Z we have that 

\\F(M)\\n = ||(L_ a (L a oFoL,))(L_ I ({a: fc })|| w 

< \\L-s\\j?( s (N )\\Ls o Fo L s (L-i({x k })\\N 

< Ci,s||i-s||if(s(No))ll^-dl^f(s(No))||a;fe||o- 



□ 



Furthermore, an order symbol matrix considered as an operator on functions, 
as in the previous lemma, extends to a bounded operator on L? (R) 

Lemma A. 5. Considering an order symbol matrix f(m,n) as an operator 

A : <S(R) S(R) 

('as in the previous lemma), A extends to a bounded operator 

A : L 2 (M) -> L 2 (M). 

Proof. This follows readily from Schur's test (for example, see Theorem 5.2 in [3]) 
but for completeness we will give the proof. Our proof is essentially the same as 
the proof in [5] , except in far less generality. 

For g e £ 2 (R), we have that g(k) := (ss^O^qr) is in £ 2 (N ). Thus 

2 



2 f£/(n,*)W) 

nGN VfcGNo / 

nGN VfcGN / 

= E f E Vmn7W\Vmn7W\\m\] 

neNo VfcGNo / 

< E f E i/(^ fc )i) IE i/(»>*)iis(*)i 

neNo VfcGNo / \fcGZ d 

<cft£ (Ed + i*- *!)-") (e 

riGN VfcGNo / \keZ d 

<c^ fx; (1+1*1)-^ f e a 

VfcGNo / \n,feeN 

VfcGNo / VfcGNo / 



(l + \n-k\)- N \g(k)f 
n~k\)- N \g{k)\ 2 
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□ 

Appendix B. Weighted L 2 bounds into weighted i°° bounds 

Now, we will prove the following lemma. It is stated in [4], but given as an 
exercise. We include the proof for completeness. 

Lemma B.l. For f £ C°° (K d ), if 

e (l + \x\ 2 f-^/ 2 L 2 (R d ) 

for all a £ 1 d , then for s' > s — d/2 and any a £ T, d 

^-{x) G (1 + |x| 2 ) (s '- |Q|)/2 L°°(R d ). 
ax a 

Proof. We may assume that \x\ > 1 on supp/, because only the asymptotic be- 
havior as | a; | — > oo matters for either property, as / is smooth. We introduce polar 
coordinates on {\x\ > 1} C K d 

x = tuj for uj £ S"* -1 and t > 1. 



9 %f(x) £ t s -^L 2 {M+ x S d -\t d - 1 dtduj). 



In polar coordinates, our assumption on / can be thus written as 

dx° 

It is clear that Xj £ tC°° (S^ 1 ), so we see that 

P-=T P-P- e t s L 2 (R+ x S d -\t d - l dtdu), 
dujk r— i oujk dxj 

3 = 1 J 

because differentiating in Xj lowers the power of t by one, but a power of t is also 
regained from the S^-, and multiplying by a function in C 00 (S' d ~ 1 ) cannot hurt the 
L 2 ^- 1 ) bounds. Similarly 

% =Ef^e t ->X>(M+x5- t '-^). 
j=i j 

Thus, we have that for any differential operator on S d ~ x , P, and any I £ No the 
assumption on / gives 

J^(P/(t,w)) £ t s - l L 2 (R+ x S d - x ,t d - l dtduj), 
We can rewrite this as being integrable on M + with values in L 2 (S d ~ v ) 

jj(Pf(t,u)) £ t s - l -^- 1) / 2 L 2 (R + ;L 2 (S d - 1 )). 

This does not depend on the exact form of P, or even its order, so we can take P 
to be very large order and elliptic, and thus for any k, we can consider original / 
as a L 2 function on R + taking values in H k (S d ^ 1 ) 

J^(/(t,w)) £ t s - l -( d - 1 V 2 L 2 (R+;H k (S d - 1 )). 

Thus, for any differential operator on S , Q of order q we have 

J^(Q/(t,w)) € t s - l - ( - d - 1 ^ 2 L 2 (R+;H k - q (S d - 1 )). 
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For large enough k, by Sobolev embedding H k q (S d 1 ) C L°°(S d 1 ), so we have 
that 



sup 



8t l 



(Qf(t,tj)) 



e t s-l-(d-l)/2 L 2 



Writing 



sup 



d 



& 



_ £S-l+p-l-{d-l)/2 



for g(t) e L 2 (R + ). Because L°° membership for smooth functions does not depend 
on its values on any compact set, g is zero for t < 2. Furthermore, 

t p ^(Q.f(tM) = J* A (t*^(o/(f, w) )) df 



so 



sup 

tGR+,u)GS <i - 1 



t»|i(5/(i,^)) 



i: 



dt' 



s-;-(d-l)/2+p-l 



< 



,2s-2i- 



-l+2p 



1.9(f) I*' 

1/2 



df 



\g(t')\ dt' 



1/2 



This is finite for 2s — 21 — d — 1 + 2p < —1, or in other words s — I — d/2 + p < 0. 
Reverting to regular coordinates, this proves the lemma. □ 

Appendix C. Boundedness on L 2 implies a bounded symbol 

In this appendix we prove that given a isotropic pseudodifferential operator of 
order close enough to zero which is bounded on L 2 necessarily has a bounded 
symbol. The proof basically involves the idea of a pseudodifferential operator acting 
on an exponent, but we could not find this particular setting in the literature. 



Lemma C.l. For 



for some e < 1/2, if Q is bounded on L 2 (R), then 



the symbol of Q, ctl(Q) — q(x,£) is bounded, i.e. 

sup \q{x,£)\ < co. 

(i,Q6l 2 



Proof. Because Q is bounded on L 2 ( 

WQfh'c. 



) there is some C > such that 

< VcwfWw 



Choose e > 0, z € R and (p(x) e C^°(K) a bump function with ip(x) € [0, 1] and 
if = 1 for \x\ < e/2 and ip = for |a;| > e. Let p zo (x) = <p( x ~ z o)- For v G [—1, 1] 
we thus have that 

\\Q(e*y 2 ^ Z0 (y))\\h m < C\\ Vzo \\h m = C||^||| 2(R) . 
Thus, we have that 

(C.l) I v (z ) := ( e -^g*g(e^>(y)),^) L2(R) < C|M|£ a(R) . 
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Letting p(x,£) = <tl(Q*Q) 

I v (zo) = J e^-^p(x,Oe i{y2 - x2)v Vzo(y)^)dyd^dx 

= J e i( - x - y ^p(x + zo, £ + 2vz )e i ^- x ^ v ip{y)ip{x) dyd^dx 

= J e^ x - y ^p(x + zo,2vx + 2vzo + r 1 )e^ v - x ^ v tp(y)tp(x) dydrjdx. 

In the last line, we made the substitution £ = 2vx + rj. Now, expanding around 
7] = 

p(x + zo, 2vx + 2vz + r]) 

= p(x + z , 2vx + 2vz ) + / p 2 {x + z , 2vx + 2vz + trj)rj dt. 

Jo 

Combining this with the above, the first term simplifies via Fourier inversion giving 

p(x + zo, 2vx + 2vzo)p(x) 2 dx. 



Now, we will show that as \zq\ — > oo the second term goes to zero. Let x € Co° be 
a bump function x(v) € [0, 1], xiv) = 1 f° r \v\ < 1 an d xiv) = f° r \v\ > 2- Using 
this we can estimate the second term where r] is small 
r t=i 

e i(*-y)v p2 ( x + ZQj 2vx + 2vz + tt]) 



t=o 



VX(v) e Vl P(y) l f( x ) dtdydr/dx 

-J J \P2(x + z ,2vx + 2vzo+t?])\\T]\\x(v)\\ ( p(y)\\ ( p(x)\ dtdydtjdx 

< C J (1 + \x + z \ + \2vx + 2vz + Hf - 1 

H\x{v)\W{y)\W{x)\ dtdydrjdx 

^'(l-e+N) 60 " 1 . 

In the last line, we used that on the support of ip(x), \x\ < e, so 

\x + z o\ > \ z o\ — \ x \ > \zq\ — e 

which gives the above inequality. Similarly, when r\ is bounded away from zero, we 
can integrate by parts, obtaining the bounds (using v G [—1,1]) 

/ ft o J^ d y (el{x ~ v)v)p2(x + Zo ' 2vx + 2vz ° + tr,) 

77(1 - x(v)y (y ~ x)2v f(yMx) dtdydrjdx 



IS, 



t=l 1 



( e *(x-y)n) p2 ( x + z 2vx + 2vz Q + try) 

=0 T 

r?(l - x(v))d*(e i{y - xfv tp(y)Mx) dtdyd V dx 
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<C"(l-£+|zo|) e °- 1 f f 1 11 ■ f 2 (7?)l l^ 3 (e % " x)2 >(y))lb(^)l dtdydndx 
J Jt-o \m 

<c"\i-e + \z \y°-\ 

Combining all of this into (|C.1|) we thus have 

(C.2) J p(x + 2d, 2vx + 2vz )tp(x) 2 dx + R(z ) < C|M|| 2(H) , 

where R(zq) is the sum of the above two terms, and by the above bounds, we know 
that R(z ) — > as \z \ — > oo. Because all derivatives of p go to zero, p is uniformly 
continuous, and thus for S > 0, we can take e in the definition of ip small enough 
so that 

/[ sup [p{x' + Zo, 2vx' + 2vzq)] — p(x + Zq, 2vx + 2vzo))ip(x) 2 dx < S. 
i'£B ( (0) 

This combined with (|C.2j) gives bounds of the form 

\p(z Q ,2vz )\ < P 

for some P > (depending on e) and for large \zq\ (and thus for all zo because p is 
certainly bounded inside of a compact set). From the above proof, it is clear that 
P does not depend on v £ [—1, 1], so we know that 

sup \p(x,£)\<P 

\x\>2\£\ 

by taking appropriate v and zq in these bounds. For q the symbol of Q, this shows 
that 

sup |g(x,C)l < V^P. 
M>2|?| 

Now, to extend this to all of R x 1, taking the Fourier transform of the operator 
Q, denoted Q and applying the above argument to Q (which is still bounded on L 2 
because the Fourier transform is an isometry L 2 (M.) —> L 2 (R)) we have that 

sup \q(£,x)\ < VP 7 . 

|x|>2|€| 

This shows that q is bounded. □ 
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